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o 

^^ , Ising Monte Carlo simulations of the random-field Ising system Feo.80Zno.20F2 are 

^vj ' presented for _ff = 10 T. The specific heat critical behavior is consistent with a ~ 

and the staggered magnetization with /? ~ 0.25 ± 0.03. 
^: 

It has recently been shown experimentally ||] and through Monte Carlo (MC) simulations 
that the three dimensional {d = 3) dilute, anisotropic antifcrromagnet Fea;Zni_a;F2 in an applied 
uniform field exhibits the equilibrium critical behavior of the random-field Ising model (RFIM) p| 
only for magnetic concentrations x > 0.75. Although there is some agreement between previous 
MC and equilibrium experimental data, particularly the staggered susceptibility critical exponent, 
^ . where 7 — 1.7 ± 0.2 for MC H and 7 — 1.58 ± 0.08 for experiment M, and the correlation length 

critical exponent, where v = 1.1 ± 0.2 for MC M and v = 0.88 ± 0.05 for experiment [Q, the 
experimental value B of the specific heat exponent, a ~ 0.00 ± 0.02, disagrees substantially with 
the simulation result 0] a = —0.5 ± 0.2. The order parameter exponent, /3 = 0.00 ± 0.05 from 
simulations 0], has yet to be measured experimentally. The MC exponents violate the Rushbrooke 

F^ \ scaling inequality 2/3 -|- 7 + a > 2, indicating possible inconsistencies. The previous MC studies 

r^ . employed finite scaling on ferromagnetic lattices of sizes L x L x L, with L < 16. Simulations 

O ' were made over many random-field configurations and equilibrium was ensured by insisting that 

, ^, . different starting configurations yielded consistent results at each temperature. In the present MC 

study, we have explored a much larger antiferromagnetic lattice and modeled our simulations closely 
after Fei_2:Zna;F2 and the thermal cycling procedures used in experiments. Our large lattice size 
^ ' prevented the averaging over many random-field configurations, though the results obtained for a 

^z^ . few different configurations yielded consistent results. We determined the staggered magnetization, 

p^ ' Ms, and Cm versus T. The resulting Cm more closely mimics the symmetric, nearly logarithmic Cm 

Q^^ , peak of the experiments Q than the nondivergent cusp found in the earlier MC studies. In addition, 

^^ ' we find (3 = 0.25 ± 0.03, in disagreement with the previous MC studies. 

^P , The body-centered-tetragonal magnetic lattice of Feo.80Zno.20F2 is modeled as two cubic sub- 

lattices delineated as one-dimensional arrays bit coded to allow for a very large lattice size, 2LxLxL 
with L — 128, corresponding to 3.4 x 10^ spins. Periodic boundary conditions arc imposed. The 
dilute MC Ising Hamiltonian is 
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^ I i7 = J2 ^ EiE-jSiSj - h'^EiSi (1) 

o 

O ■ where Si — ±2 and e^ = 1 if site i is occupied and zero if not, J2 = 3.17 K and h = 8.73 K. This 

approximates the experimental system which has one dominant Heisenberg exchange interaction, 
J2, and a large single-ion anisotropy. The very small, frustrating interactions, Ji and J3, are not 
included. The MC value for J2 is 0.60 times that of FeF2 so that T'c(O) corresponds roughly to 

C^ ' that of the experimental system. The value of h is chosen to have roughly the same effect as a 

field iJ = 10 T in the real system, as determined by the fields at which low temperature spin flips 
occur. This field also yields a shift Tc(0) — Tc{H) consistent with iJ = 10 T in the real system |^]. 
Henceforth we refer to the applied uniform field as iJ = 10 T. 

For H — 0, the sample is cooled to low temperatures in steps of 0.01 K while magnetic sites are 
randomly visited an average of N times, where 2000 < N < 5000, and flipped with a probability 
given by the metropolis algorithm. Similarly, the sample is heated from an ordered state to above 
the transition. No hysteresis is observed at iJ = using these two procedures. For H > 0, 
two thermal cycling procedures are used to mimic those of the experiments. Each site is visited 
5000 < N < 10000 times per 0.01 K temperature step. For zero-field coohng (ZFC), the sample is 



started at if = in a well ordered state (in experiments this state is prepared by cooling in iJ = 0), 
the field is subsequently applied, and the sample is heated through Tc{H). In field-cooling (FC), the 
sample is cooled with H applied. We observe no significant hysteresis with these two procedures, 
consistent with the experiments. Ms and Cm are calculated for each temperature. The simulations 
are performed on Linux/Gnu computers with CPU speeds between 350 to 450 MHz. Individual runs 
with L = 128 at 5000 MC steps per spin and covering lOK in T typically take one week of CPU 
time. We have not seen any significant dependence on N in the range studied. 

Figure 1 shows C™ versus T for L = 128 and N = 5000 for i7 = and for ZFC and FC at 
H = 10 T. Each point represents a point-by-point derivative of the energy which is itself averaged 
over 0.1 K intervals. The behavior is very similar to that observed in birefringence and pulsed heat 
experiments Q. Namely, the H = behavior exhibits the asymmetric cusp (a < 0) characteristic 
of the random-exchange Ising model and, for 7J > 0, a much more symmetric peak near Tc{H) 
consistent with the symmetric, nearly logarithmic divergence. Previous MC studies employing finite 
size scaling did not yield this behavior. 

Figure 2 shows the ZFC M^ versus T - Tc ioi L ^ 128 and N = 5000 for iJ = and i/ = 10 T. 
Fits of the data for various ranges of |i| within 10""^ < \t\ < 10^^, yield the critical order-parameter 
exponent /3 = 0.35 ± 0.04 for 7J = (REIM) and /? = 0.25 ± 0.03 at iJ = 10 T (RFIM). The TdH) 
values used in the fits were constrained by the fits to Cm{H) 

To demonstrate that the results of the simulations are close to equilibrium, we did a simulation 
just below the transition for L — 128, H = 10 T, and T — 60.5. One lattice was started fully 
ordered, another fully ordered but with the sublattices reversed, and one randomly ordered. All 
three lattices converged to the same energy and Mg within expected fluctuations after 5000 steps 
per spin. 

Using P = 0.25 ± 0.03 from MC and a = 0.02 ± 0.02 and 7 = 1.58 ± 0.08 from experiments, 
we obtain 2/3 + 7 + a = 2.08 ± 0.14, which is consistent with the Rushbrooke scaling relation as 
an equality. Further studies will be needed to understand why the MC simulations with large an- 
tiferromagnetic lattices and thermal cycling yield different results than the earlier ferromagnetic 
MC simulations with finite scaling analyses . For simulations of RFIM critical behavior employ- 
ing dilute antiferromagnets, the magnetic concentration should be kept well above the percolation 
threshold concentration for vacancies. 
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FIG. 1. Cm vs. T for ii" = and for ZFC (triangles) and FC (inverted triangles) at 10 T with L = 128 
and X = 0.80. Note the symmetric peak for H = IQ T. The curves are fits to the data with a — —0.09 and 
0.00 for the H = and H = lOT cases, respectively. 
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FIG. 2. Ms vs. T - Tc for H = and ZFC at H = 10 T with L = 128 and x = 0.80. The curves are fits 
to the data with l3 — 0.35 and 0.25 for the H = and H — lOT cases, respectively. 



